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ABSTRACT 

This study considered the feasibility of obtaining 
a self-adapting pitch rate control system of a hypothetical 
airframe, utilizing linear compensation of the invariant, 
i.e., normal operating state, parameters of the aircraft. 

The relationship between the flight dynamics and 
aircraft parameters was established, and each was reduced 
to an invariant and variant component. A method of repre- 
senting the system with, and without parameter variation 
was then developed. Linear compensation of the normal op- 
erating state was developed to provide nearly flat, minimum 
phase shift, frequency response; as the primary performance 
criterion. 
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The compensated and uncompensated aircraft was 
analyzed in widely varying operating conditions. Theoret- 
ical frequency response and analog computer transient 
response simulation, yielded the following results: 

The desired frequency response in all flight con- 
ditions can almost be obtained, since the compensated 
system Introduced varying attenuation with a maximum 
of 35^. Six of the seven conditions produced acceptable 
transient response. 

To eliminate the attenuation, a variable gain 
amplifier, sensitive to the conditions is proposed. Us- 
ing linear compensation the complete system, including 
the model, amplifier, and pitch rate system, acceptably 
reproduced the desirable model output, through most of 
the various flight conditions. This satisfies the 
basic requirements of a self-adaptive system. 

Thesis Supervisor: Y.T. Li. 

Title: Associate Professor of Aeronautics and 

Astronautics . 
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OBJECT 



The object of this study is to propose and 
investigate a self-adaptive control system by 
means of compensation of normal operating state 
aircraft parameters, using linear techniques. 
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CHAPTER 1 



INTRODUCTION 



Modern Control System Requirements 

With the advent of high performance aircraft# 
satellites, rocket aircraft, guided missiles, ICBM's, 
and eventually, space ships, there has been an acceler- 
ated requirement for replacement of the human pilot with 
automatic flight control systems. This replacement is 
not necessarily complete, and in many cases is used only 
to augment the human characteristics in order to add damp 
ing, Increase accuracies, and in extreme cases, to provide 
control capability. This class of automatic control 
system has been further instrumented, and is more popu- 
larly known as a self-adapting control, or servo, system. 

The human servo system is an excellent example of 
self-adaptation. However, in the present day supersonic 
vehicles, reaction time is too slow, and muscular coordln 
atlon is not sufficient to cope v#ith the various extremes 
of velocities, temperatures, forces and other modifying 
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inriuences; since, modern aircraft have flight capa- 
bility over extreme ranges of altitude, air speed and 
Mach niimber. This capability is mostly attainable due 
to propulsive thrusts which approach maximum aircraft 
gross weight. Automatic control systems and stability, 
augmentation devices have thus been installed in aircraft 
to provide satisfactory stability and control throughout 
the various flight regimes. The stability and control 
characteristics of all aircraft vary as functions of 
equivalent air speed, Mach number, altitude, aircraft 

gross weight, center of gravity position and external 

( 1 ) 

configuration. 

System Representation and Classification . 

Self-adaptation may often be represented by an 
ordinary negative feedback control loop which in Itself, 
is the first of the basic servo concepts. However, the 
quality of this adaptation may, then be used as a measure- 
ment of performance capability. 

Some basic classifications of the adaptive concept; 
The first, and most easily understood is that of negative 
feedback. If there 1s a high gain in the forward loop, 
assuming non unity feedback then the closed loop transfer 
function then becomes approximately that of the negative 
feedback performance equation. The forward loop has no 
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effect essentially, on the entire closed loop if the 
combined gains are much greater than one. Thus, if 
the forward loop is the controlled member, then the 
system operation will be completely independent of 
the controlled member, e.g. , an aircraft, which is 
an ideal situation in self-adaptation. However, 
actual implementation with mechanical or electrical 
devices are not as easily found as the transfer 
equation would have us believe. 

There are many ways of re-writing the block 
diagram representing the above concept in order to 

overcome the basic difficulties with normal feedback. 

( 2 , 3 ) 

The method that is Introduced by a number of authors 
is that of the conditional feedback loop based on a 
performance model, that is designed such that the over- 
all system approximates the characteristics of the 
idealized model, over wide ranges of operation. 

( 2 ) 

Sperry Gyroscope Company have further described 
systems such as conditional feedback with a pre-filter, 
which will be further discussed in the later chapters 
of this report. 

Compensation and Effects . 

Automatic control systems are generally multiple 
closed loop with feedback paths utilizing variously: 
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displacement, velocity, and acceleration of the aircraft 
as control quantities. Therefore, the static and dynamic 
response characteristics have a great effect and cannot 
be ignored by using only high 'gain feedback, since satura- 
tion and power obtainable from the physical system often 
falls far short of the theoretical desires. It becomes 
obvious, then, that some compensation must be made for 
the aforementioned effects. 

This compensation may take many forms as predi- 

( 3 ) 

cated by V^Tiitaker, Yamron, and Kezer;, namely, changing 
system sensitivities as calibrated, programmed functions 
of indicated airspeed, Mach number and other flight quanti 
ties. The adjustments to the system sensitivities are 
then made such that they are open loop, from the stand- 
point of the effect of the adjustment on system perform- 
ance, Fhst action adapting control systems are now be- 
coming of greater importance, as stated before, and will 
become even more so with the advent of space flight. 

There will be colossal changes in environmental parameters 
as well as the basic system parameters themselves. Large 
changes in mass with attendant changes of inertia, aero- 
dynamic heating due to high velocities in the earth’s 
atmosphere, and changes of the density of the mediiam in 
which the controlled vehicle is operating, are a few 
examples of the changes that must be anticipated and 
corrected for by the vehicles control system. Recently 
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it was announced that the problem of re-entry of a 

manned satellite would be undertaken. Let us assume 

that the satellite would re-enter at an altitude of 

100 kilometers, and utilize a bounce-glide technique 

of re-entry, i.e., decelerate gradually in the upper 

layers of the sensible atmosphere until a tolerance of 

temperature is obtained, then bouncing back to the next 

higher level and then descending again, in gradual 

sweeps, for re-entry. This requires that a theoretical 

aircraft, weighing 20 tons with a wing loading of 30 

kilograms per square meter, undergoes skin temperature 

variations of from 800 degrees G, to 400 degrees G, 

to 1400 degrees G, to 0 degrees G, while undergoing a 

velocity change from 8,5 kilometers per second to zero 

(4) 

during a 2 -hour flight. 

It is seen, then, that it is necessary either to 

change certain system parameters; or to modify command 

signals, as the environment changes, in order to maintain 

satisfactory dynamics. It is further necessary, since 

flight testing in outer regions of space vehicles is not 

physically realizable, to be able to design a theoretical 

system that can be implemented at a later date that is 

intelligent enough to know when dynamic performance has 

changed and to utilize this knowledge to a dap t itself, 

so that its final performance will be within a satisfact- 

(3) 

ory range at all times. 
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The present trend to supplement, or even eliminate, 
a human pilot, has led to severe requirements for control 
performance, since the capability of many applications re- 
quires control through violent maneuvers and yet must still 

( 10 ) 

be sensitive to small input command signals. .In many 
cases a control system may be extremely complex; however, 
when dealing with a specific vehicle such as an aircraft, 
we are limited by size, weight and structural consider- 
ations. The result, then, is often a servo with a limited 
output rate. If we are able to make this assumption, the 

controller may be considered as a simple linear system 

( 14 ) 

with a limit on the rate of control. But, it is this 
very limit on the control surface rate that is a major 
difficulty in a system requirement for stability for both 
large and small commands. Compensation and adjustment 
of the system parameters to give the required stable res- 
ponse generally results in undesirable slow response to 
small signals. There have been many theories on compen- 
sation for such undesirable effects, more specifically: 
to design system parameters as functions of error and 
error rates. This has led to development of non-linear 

elements to improve the performance of both linear and 

( 5 ) 

non-linear devices.. 

There have been methods proposed to also accomp- 
lish this, in which a control system tends to adjust 
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itself to a pre -determined operating condition. The 
system Imposes known continuous variations on system 
inputs to determine system parameters, to maintain out- 
put at a maximum. 

Series compensation and unity feedback with a 

controlled component in which the series compensation 

is varied by an external computer is yet another propo- 

( 6 ) 

sal. 

Summary 

Thus, adaptive mechanisms may be cataloged into 
two general types, the first of which is the simple high 
gain system using linear negative feedback. The second, 
which is of primary Importance in the recent history of 
self-adaptation, is where the alteration of the control- 
ler as a function of system response, is involved in a 
non-linear system. This type leads ultimately to an op- 
portunity to design and manufacture a controller that 
is optimum Insofar as the continuous detemination of 
the static and dynamic characteristics of the controlled 
members is concerned. 

To correct the difficulties in stability already 

mentioned, it has been s uggested that uniform stability 

may be obtained by compensating for changes in aircraft 

( 12 ) 

stability directly. This is accomplished by adjusting 
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control parameters as functions of response in order 

( 8 ) 

to maintain a unlfom response. The basic instabilities 
generally result from high feedback gain, due to the 
dynamic effects of the servo actuators. Todays manufact- 
uring techniques, however, have been able to synthesize 
controllers with better rates of response, and thus 
higher gains are obtainable, leading to more practical 
value in the theory of negative feedback. 

By definition, optimum control of the aircraft 
results if the control system provides close control of 
the transient and steady state responses of the aircraft 
to command and disturbance inputs, so as to use the 
capabilities of the aircraft most effectively for ful- 
filling the specifications of its flight mission. The 
characteristics involved being response time, damping, 

dynamic and static errors, and control of interference 

(2} 

effects. This very definition allows great flexibility 
in utilizing the aircrafts capability, without forcing 
an impractical response from the physical vehicle. There 
are as many design criteria as there are designers, how- 
ever, and in selecting such a criteria it must be remem- 
bered that the desired performance itself may vary from 
one phase of a flight to another. 
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CHAPTER 2 



SCOPE 



The purpose of this study is to Investigate the 
feasibility of obtaining a self-adaptive system by em- 
ploying a linear technique, herein defined as invariant 
parameter compensation. The basic philosophy underlying 
this technique is the estabilishment of a state in which 
the airframe parameters, and hence the airframe dynamics 
are defined to be invariant. Phrther, having established 
the airframe in this state, it is then assumed that suit- 
able compensation can be devised to produce an essentially 
flat response within specified tolerances, of the air- 
frame controlled system over a wide range of frequencies. 

The design tolerance criteria is as follows; 

a. Amplitude ratio l:jtl^ variation 

b. Phase shift of less than ±10® 

all within an excitation frequency range of 0-50 
radians/ second. 

The value of this technique lies in the capability 
of establishing an effective natural frequency very much 
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greater than any frequency to which the system will 
normally be subjected, and an effective damping ratio 
which is dynamically acceptable. Through this technique, 
it is believed that the normal range of parameter varia- 
tions to which the aircraft will be subjected will not 
change the performance function to any great extent, 
over the range of interest of input frequencies. It 
is further desirable that the aircraft perform in a 
manner as specified by the performance of a reference 
model system which precedes the controlled loop, and 
whose output will be the input to the cbntrolled loop. 

The performance function of the model will be developed 
in later chapters. 

If the aircraft, as compensated, does reproduce 
as its output, the output of the reference model system 
for all normal variations of parameters, and for all 
frequencies within a range of interest; then the aircraft 
control system will truly be self-adaptive. 

The scope of this paper will be limited to the 
pitch rate perfomance of a hypothetical aircraft, whose 
performance function and dynamics are; 



Performance Function Diagram; 

AIRFRAME 





Qy. 


Phi. 


S. 
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Whe re : 



PI^ — Reference model system performance function, 

PPs^= Servo and Actuator performance function, 

Za= first order zero; an airframe dynamic, 

Wj^= Airframe natural frequency; an airframe 

dynamic, 

:f = Airframe damping ratio, an airframe 

dynamic, 

0 

0 in =• Commanded pitch rate. 

©m= Model commanded pitch rate. 

S E=Elevator deflection .angle. 

a 

0o =Airframe pitch rate response. 

In order to employ this method of obtaining self- 
adaptation, it will be necessary to establish the relation- 
ship between the above airframe dynamics and their correspond- 
ing parameters, Phrther, since compensation is to be performed 
on only the Invariant parameters, it is necessary that a 
method be devised wherein the invariant parameters may be 
handled separately from the variant components. 

This development will lead to representation of 
parameter variations by synthesis of a fictitious system 
containing both the invariant and variant components. The 
variant components represented in this fictitious system 
are the inner loops, which are active only when the air- 
craft is not in its invariant state. 
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Although it is anticipated that this technique 
can be utilized by any aircraft configuration, a hypoth- 
etical airframe was employed, since by this method the 
contributions of the various parameters to the aircraft 
dynamics and the effect of variations of these parameters 
on the dynamics could be more adequately demonstrated. 

The invariant system parameters are those which 
exist when the subject aircraft is in the state of alti- 
tude, airspeed, . and loading specified as the invariant 
state of operation. The specification of this state appears 
to be. somewhat arbitrary, although it would seem logical 
to designate the ’'normal" operating state of the aircraft 
as its invariant state, Ibr aiframes not possessing such 
a state, selection of the invariant parameters becomes 
a matter of additional study and specifications of criteria 
to establish the "best" set of invariant parameters for 
such types of vehicles as ballistic missiles, may well 
prove to be of sufficient magnitude to justify further 
work in this field. 
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CHAPTER 3 



THE HYPOTHETICAL AIRFRAME 



Definition of Parameters 

As indicated in Chapter 2, the effectiveness of 
linear compensation is to be tested on a hypothetical air- 
frame. In this chapter the airframe pitch rate equation, 
dynamics and parameters are established. In addition, the 
invariant state will also be defined. 

For the purpose of this development, the following 
table of parameters is applicable: 

Dimensional 

•=r air density 
= air density at sea level 
— airframe airspeed 

= a characteristic surface area of airframe 
= a characteristic chord length of airframe 
= moment coefficient 
= coefficient of lift 

= derivative coefficient = dCl ; i =m,L ; j = 



f 

fsL 

V 

s 

c 

C»y( 

Ct 
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m 









Ce = 
o< = 
G = 
Q = 



iyy= 




m 



control surface deflection 
airframe angle of attack 
airframe pitch angle 
airframe pitch rate; l,e. d0 



longitudinal axis moment of inertia =k 
longitudinal axis radius of gyration 
mass of airframe 



Non-dimens lonal 



(T 






c 



Aircraft Performance rhnction 

The short period performance function relating 
elevator angle (Of) input and aircraft pitch rate (0©) 
output may be represented in terms of dynamics as follows: 

4 k[^ 

(Note; The negative relationship is introduced 

( 11 ) 

by sign conventions for Oe f 6© • should be apparent 

that the same sign convention will apply for ( Qt , thus 
making a positive relationship for 0^° y • ) 
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It Is assumed that each of the dynamics my be repre- 
sented as the sum of a variant and an invariant component, 
as follows ; 

Za = Zao + A Za ; K = Ko 4- ^ K ; /Wn = (TWn)Q+A (/Wn) 

(wg)^ + a(w2) 



The subscript zero signifies the parameter is in- 
variant. The selection of & j*V/n as basic parameters 
was done to eliminate cross product terms in later develop- 
ment. 



Dynamics in terms of airframe /flight path parameters 



W' 



I = where; 



A — C»V>< 



where; 



i fm - 

2a = A V(T 

^s/< *A^(r 

2 

^ = Al V where; 



— C/»o 



j A2 - 

Ku ^ J 



-fc Ku 



where ; 1 ] 

-'i CyriS^ ; A5 = C»,cJ 









Aq = 



The following parameters are assumed to be constant: 
2 



S, C, Ky» > 



C„,e( » C^o< » i 



& C* 



» ~ yy-jV, » -toj » -/nv, » ® ACtUally, 

the derivative coefficients vary somewhat with Mach number. 
However, the variations of these quantities are not signifi- 
cant when compared to variations of other parameters. Conse- 
quently, the results obtained will not be too great in 
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error as a result of this simplifying assumption. 



Just as the airframe dynamics are considered to 
be composed of a variant and an invariant component, so 
too, may the airframe/f light path parameters be considered. 
The following then define the dynamic components in terras 
of the corresponding parameter components: 



(a) ^ (jJ^ 

^ ^ A (7T 1^ 

A (jui = ■ —I ^ + 

M2[(<^ f A<r)Vv.tAV)-(c;v;f] - { u>.)l (z ^ ^X)] 



(c) Ze. = 



(d) 









= i 

(<5r^A(T) 




K - K> ^ ^ K 

= .^7 ^^(TqXo 1-AnAr Vj" (T: 

f R ^ f C") - (Vo J 

- /<o ^ 
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specification of the Invariant State 



The hypothetical airframe is assumed to be at its 
invariant state, i.e. normal operating state, at an alti- 
tude of 35,000 ft., with an airspeed of 1000 ft. /sec. and 
a mass of 497 slugs. This condition produces the following 
values, defined to be the invariant values, of the parameters 
^ _ .835 X 10^ slug/ft3 (T^ = *35 

Vo lOOOft/sec. ^ 4 ^ = 79.1 (see below) 

Mo = 497 slugs 

In addition, the parameters previously defined as 
cons tants have been assigned the following values: 



Parameters Coefficients 



Ky 


81 ft^ 


Al 


— 


-0.00275 


Ky 


1.046 


^2 


— 


0.0485 


s 


300 ft^ 


^3 


= 


-0.475 


c 


8.8 ft 


A4 




-0. 486 




-4.5 


As 




-2.0 




-.45 


As 


= 


-0.6 


Qwtn 


-3.0 


^7 




0.00306 





-7.0 


^8. 


= 


-0.30 


II 


-0.2 








Q»>^' 


-0.5 
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utilizing the foregoing parameters, the invariant 
dynamics are determined to be: 

(V<n>o= 

2(fVVn)o= 2.1 

(Za)o = 1-07 

(K)^ = 13.55 

Once the dynamics have been defined, the problem 
is reduced to one of compensation; to obtain acceptable 
aircraft performance. 
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CHAPTER 4 



COMPENSATION AND EPPECT 
OF PARAMETER VARIATIONS 



Method 



As previously indicated, this study is concerned with 
compensating the air frame for one state of operation and 
then testing the effectiveness of this system when the vehicle 
is in states other than that in which it is compensated. 

It was found, that handling of variational effects could be 
greatly facilitated by employing the analog representation 
of variant and Invariant parameters shown below. 



(a) Phnctional Diagram and Exact Equation of 

Of 

Diagram of elevator to aircraft 




Define KG^^'q such that: 

(P^+ ^JWn p + Vl/^) qo = K. (p -I- Za)<$£ 

(1) [p2+ 2([fo;^\+ A/4)p+(wg)o+A(W^}] qo = 

(Ko+AK) (p-4- Zao+AZa) 

¥;here the subscript o refers to the invariant 

parameters. 
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(b) Analog Representation of Equation (1) 




(2) q^= KoGo (1+ >4 -(KoGoHp) q^ 

where KoGo = Ko (p Zao) 

p®i (Wglo 



select: G^^ , anct Hp such that substitution 

into (2) yields (1) 

let: G^ = ^ Za ; ( P f Zao +a Za ^ ; 

(p+Zao) Ko (pfZao) - 

Hp=.2 A (fVvn)p- A(W§) 

Ko (p Zao) 

then; (p^-*- 2 f-/£4,|p + qQ=Ko(p+Zao) 

fl Za ^AK {pi- ZaofAZa)) 

\ (pi^c\o) Ko (p+Zao) j 

- {2 A(fWn) p + /i(w2)} 

or:|p2 +• 2[( ■fVVn)o ^/a(/V/n)] p f (W^) of A( V«§) j q^ = 
(KofAK) (p f Zaot^Za) 

which is identically equal to equation (1) 
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Hence, to compensate the invariant system, the 
system can be treated as though the inner loops are opened. 

Compensation of the Invariant System 

To compensate the airframe in its invariant state, 

i.' 

a simplified system was employed. The servomotor and ele- 
vator actuator dynamics are reduced and the dynamics of 
the rate gyro in the feedback loop were ignored. There 
is no loss of effectiveness caused by making these assump- 
tions. Since one of the objects is to compensate entirely 
for the servo and actuator dynamics, and the dynamics of 
rate gyros currently employed in practice are such that 
they are of small consequence when compared to the dynamics 
of the other eleraepts in the system. 



A. Simplified System: Servo & Airframe 

Actuator 




Open Loop Performance Equation 



It is desired to increase the damping ratio to a 
value of approximately .7 since this^ damping ratio produces 
desirable dynamic performance. It is further desired to 
increase to a value much greater than (W^)©, since such 
an increase will nullify, to a large extent, the effect 
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of parameter variations which occur when the airframe 
is not in its normal operating state. Rirthermore, a high 
effective natural frequency of the second order system 
ensures a "flat" response of this component over a wide 



frequency band. 



/ 



let: i \ 

“ yib t ; or the equivalent thereof. 

then : 

^ f/> <7c)^»u 

or: ^ 

f/>*- ^ ZfuJ^ V- (. P-h ^r>y 

where: ^ (uj> I i- j Z f* uuJi ^ z(foJr,)o P K 

f'=. ^('SuJ^yX ^ ZCo 

<^^(<^n.)e ■/■ fCaK^h z / ^\ / 

Choose values of/<;^ ^ such that X>^^j[and T ^ »7 



B. Determination of KcGc: 

^c) 



let: Cr 






'c 7c -V 

Then the system is effectively reduced to: 

1 



or the equivalent thereof. 



t 








J 



-f- z fcu^'y ^ 
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where : 



A frequency response analysis was made on the re- 
sultant system employing a damping ratio of .7, an effect- 
ive natural frequency of 20 rad/sec, a loop gain of 25 
and a value of Zc = 8 rad/sec. The open and closed loop 
frequency responses are shown In the chart below; 

From the frequency plots ; 



w j 


Open Loop 


Closed 


Loop 




AR (f, 


M 


N 


.2 


1 2 


.965 


0 


.5 


1 3 


. 965 


0 


1.26 


1 6 


.965 


0 


3. 10 


1 10 
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Effect of alrcraft/fllght path parameter variations 
on compensated system 



Equivalent functional diagram; 




As indicated on the first page of this Chapter, 
j^Sz. Hp are fictitious functions usea to represent 
the effects of parameter variations. The existence of such 
variations is equivalent to closing the loops containing 
these functions. With these loops "closed" the following 
analysis of the system is made. 

Performance Function: 



Unity feedback loop open; 

- /C Kc. (vo ^ 

The frequency response of the first bracketed 
term is the open loop frequency response of the compen- 
sated Invariant system. The frequency response of the 
second term is determined by substituting the functions 
for Ggi & Gji previously derived, and the resultant func- 
tion becomes; 

whe re ; ~ ( Aa t ^ K) ( ^ iaj 

-i ' 

The frequency response of the third term is 
similarly determined and the resultant function of 
becomes: — 



X + 



A 



'<2 ('^nz'P f 1) 






Kx = 



/A u)^ 

^oto 



where : 
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Thus, to produce the open loop frequency response 
of the controlled system, it is necessary to determine: 
go (jw) = AR^/Vt- 

where: ARj = log"^ { log )+ AR(i + ) - 

log AR(n-KoG>Hp)') 

'^(KoKcGdG*) ^ ( 1 f G^^-H Gfi, } -^d+KoG’Hp} 

Using this method of variations representation; 
frequency analysis plotting and transient response, by 
analog simulation, were made for the compensated airframe 
in seven different flight conditions. 

The flight conditions tested were chosen as a re- 
sult of the desire to Impose the most severe tests of the 
effectiveness of this technique of compensation. Limits 
of altitude, airspeed, and loading were determined from 
current literature, and from operational Naval aviators, 
in an effort to maintain the conditions within practical 
limits . 

Definition of flight Conditions 

Condition A represents the aircraft operating in 
its Invariant (normal) state. Condition B represents the 
aircraft operating at high altitude, at high speed, and 
with payload expended. Condition C shows the aircraft mak- 
ing an approach for landing, with payload expended. Condition 
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D simulates the aircraft at a "medium" altitude and at 
reduced speed, while fully loaded. Condition E represents 
the aircraft at high altitude, reduced speed, and payload 
expended. The P condition shows the vehicle at the same 
altitude and speed conaitions as E, but fully loaded. Fin- 
ally, condition G simulates the aircraft in takeoff condi- 
tion, fully loaded. 

The resultant transient and frequency responses of 
the aircraft in the seven test states, both in a compensated 
and uncompensated control state, for comparison purposes, 
are included as Appendices A and C. 
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CHAPTER 5 



MODEL DESIGN, VARIABLE GAIN AMPLIflER 
AND OVERALL SYSTEM PERFORIvIANCE 



Model Design 

As indicated in Chapter 2, the system design ob- 
jective is to cause the system output response to approx 
imate as closely as possible the response of the model 
to various inputs. The model is the physical embodiment 
of the design specifications of a system, and a change 
in the model is exactly equivalent to a change in system 
specifications. The first consideration, then, is to 
decide upon the system specifications, which must be 
compatible with performance capabilities of the aircraft 

The model can then be designed Independently of the de- 

(3) 

sign of the control system. 

Since the theoretical air frame was chosen with 
no specific mission in mind, the model reference system 
for this study was arbitrarily selected to be a second 
order system characterized by an undamped natural fre- 
quency of 3.5 rad/sec and a damping ratio of .7. 
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Variable Gain Amplifier 



It is shown in Appendix A that the compensated con- 
trolled system response has little phase shift from the 
controlled system input, for the different conditions 
tested over a wide range of frequencies. However, the 
system does attenuate the input to varying degrees for 
the seven test condition. 



output response duplicates as closely as possible the 
model response, be fulfilled, it is required that an amp- 
lifier be inserted between the model and the controlled 
system. However, since the system attenuates its input 
by varying amounts, depending upon the flight condition, 
it is necessary that the amplifier be a variable gain de- 
vice sensitive to the flight condition. The performance 
criterion of this device is such that: 

- 1 
1m 

Employing the following development; an acceptable 
performance, within the above criterion, of the variable 
gain amplifier is obtained: 



Where o subscript refers to the invariant state of 
the parameter. The x subscript refers to the general state 



In order that the requirement, that the system 



Define; go = ; gx= 




- gain of variable 
amplifier 
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of the parameter, and and^ are as defined in Chapter 

3. 



Required; Ky^r 



25 gx/go 



or: K 



var = 



1+25 gx/go_ 

1+25 gx/ go 
25 gx/go 



= 1 



As indicated above, the use of a variable gain 
amplifier having this performance enables the controlled 
system oqtput response to duplicate the model output res- 
ponse. 



Overall System Performance 

V»ith the reference system model and variable gain 
device established, a frequency response plot and an analog 
computer transient response were made for the overall sys- 
tem. The results of these tests are included herein as 
Appendix D. 

It was found that for all frequencies of interest 
the overall system, i.e., the model, variable gain amplifier 
and closed loop pitch rate system; can be represented by 
the following functional diagram: 
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Where Katt represents the controlled system atten- 
uation of the input to the controlled system and Kvar 

Katt = 1 

This representation was found to be valid regardless 
of the condition in which the aircraft was operating. 
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CHAPTEH 6 



CONCLUSIONS AND RECOMivEN DAT IONS 



Conclusions 



Using only the linear technique as described here- 
tofore, the pitch rate control system can be compensated 
to the extent that the aircraft response follows the ref- 
erence system model output with a phase shift of no more 
than ± 10*^ for all practicable input rates. This method 
of compensation is not completely self-adaptive, since 
the compensated system attenuates its input (the refer- 
ence system model output) amplitude in varying degrees 
up to 25^. Unless this attenuation is acceptable, a var- 
iable gain amplifier should be employed to eliminate this 
feature of the compensated system. Since the variable gain 
amplifier is programmed so that it is sensitive to the 
flight condition of the aircraft, it is an overall non- 
linear device from state to state. But, with the air- 
craft in a given state, the amplifier performs as a 
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linear element. Utilizing the combination of linear 

compensation and a programmed variable gain amplifier, 

a pitch rate control system having essentially a "flat'* 

0 

frequency response, i.e.,J:10 phase shift and an ampli- 
tude ratio of 1 ^ 1%, is feasible. 

As indicated in the scope, the basic objective 
of this method of compensation is the attainment of a 
flat response. In other words, the system was designed 
with frequency response characteristics as the principal 
consideration, fhis system when so considered became 
self adaptive, to a great degree of accuracy. 

Of secondary consideration in this design, but 
of great significance in actual practice, are the transient 
characteristics of the compensated system. Ideally, the 
transients of the compensated system would be of small 
significance relative to the reference system model tran- 
sients, or at worst, no greater than the model transients. 
Under this ideal arrangement, the response time of the 
overall system used in this paper is the response time 
of the second order model, i.e., .85 seconds to come 

within 5% of correspondence. The analog computer tests 
of the overall system. Appendix D, shows that this sys- 
tem does approach the ideal transient performance. How- 
ever, the computer results are misleading as a result of 
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mechanization of the problem, in that the computer 
analysis was made for the general case and conditionso 
The general case was then extrapolated by changing of 
dial settings only, v/hich did not allow the required 
flexibility for proper solution to flight condition P» 

The step response equations representing the 
air frame dynamics in its various operating states as 
tabulated in Appendix C, show that with the aircraft 
in flight conditions A, B, and C, the transients are 
of no consequence relative to those of the modelo Con- 
ditions D, E, and G yield transients, which are of the 
same order of magnitude as those of the modelo With the 
aircraft in condition P the response time of the control 
system is noticeably greater than that of the reference 
modelo Aoaitional study of performance aquations revealed 
that the minimum value of the quantity required to 
produce transients which are iaeal, within the previous 
definition, is 2o3» All of the flight conditions tested, 
except P, have values of this quantity which are greater 
than this minimumo 

Exclusive of the above limitation, the method of 
invariant (normal state) parameter compensation is effect- 
ive in obtaining self adaptation in the pitch rate control 
system. It is reasonable to presume that this same technique 
will be valid and useful in obtaining the desirable level 
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of self-adaptability, in any control system 



Recommenaations 

It is recommended that; 

lo The system as developed in this study be mechan- 
ized and tested furthero In the mechanization of this 
system, it is envisioned that a rate gyro, such as the 
MIT 10^ rate gyro, with a static sensitivity of 1 mv/mr/sec 
be employedo Such a device ensures the attainment of gain 
requirements without saturation problemso To attain certain 
of the performance functions employed in the development 
in Chapter 4, it may be necessary to utilize high gain 
devices fed back upon themselves through devices having 
transfer functions which are the reciprocals of those 
desiredo 

2o As an intermediate testing and analysis procedure, 
it is evident that another type of analog simulation of 
the linear adaptive system should be made, using basic 
mechanisms and associated electronicso The characteristics 
of these aevices should closely approximate those specified 
in recommendation 

Pbrther, some of the amplification factors and 
controls may be simulated by simple hydraulic valves and 
connections . 

This "breaaboard” type of study would satisfy a 
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basic design goal; that of using available components that 
are already capable of design and/or manufactureo These 
devices, once proven, can be extrapolated and modified 
as designs for the specific devices to be used in the 
control system desiredo 

3o Additional study be made on this technique of 
compensation with the objective of obtaining a compromise 
system which yields an acceptable frequency response and 
an acceptable transient performance, ioeo, near ideal, 
for all operating states of the aircrafto 
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APPENDIX A 



Frequency Response of the Pitch Rate Systenio 

lo Definition of the Basic Plight Conditions o 

2o Polar Diagrams of the Compensated Systemo 

3o Polar Diagrams of the Uncompensated System^ 
Unity Ffeedback onlyo 



37 



lo The Basic flight Conditions Defined? 



Condition Altitude (ft) Airspeed (fps) Mass ( slugs ) 



A 


35,000 


1000 


497 


B 


70,000 


1000 


248 o 5 


C 


Sea Level 


500 


248o5 


D 


35,000 


500 


497 


E 


70,000 


750 


248 0 5 


F 


70,000 


750 


497 


a 


Sea Level 


250 


497 




2. (ARSR) 


‘lout 

Imodel 





38 



I 



I 

n 

I 

I 




I 

I 

I 


















I 





f 




J 





APPENDIX B»1 



The GAP/R Analog Coraputeri General 
Information & Techniques » 
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General Information 



The G.A. Philbrick Researches Incorporatea, Elec- 
tronic, high speed analog computer was used to verify 
the transient response characteristics of the uncompensated 
aircraft pitch rate system with unity feedback, the com- 
pensated system, and the idealized model transient response » 

It was decided to test the pitch rate system, un- 
compensated and with unity feedback only, since this is 
an obvious first approximation to self adaptability, and 
then compare this with the response of the final compen- 
sated unity feedback systemo 

The Computer used is one of the earlier models of 
the GAP/R unito The idea involvea is that of a combination 
of "black boxes" or modules, that perform specific functions 
within themselves, such as multiplication, addition, inte- 
gration, differentiation, squaring, and simulation of non- 
linear functlonso Rirther, these capabilities may be ex- 
tended to include complete first and second order servo 
systems. The modules were designed to operate in "fast" 
time, so that results may be seen on any standard oscil- 
loscope trace, from which permanent photographic records 
can be made. In the modules, herein after referred to as 
boxes a that are concerned with time functions, there are 
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characteristic time constants (T^) settings of 4 and 
o4 millisecondSo Modifications were made at MIT that 
allowed "real" time operations, by changing the value 
of the condensers in the box that is concerned, so that 
4ms, o4ms and 1 second, are available for Tq settingso 

The boxes actually used in the computer Simula” 
tion are described more fully in the latter section of 
this appendixo 

The nominal voltage range for all variables is 
100 volts, that is, minus 50 to plus 50o The manufact- 
urer specifies that electronic variations affect the 
sensitivities by less than 1^» Each box is made up of 
one or more DC minimum drift amplifiers that have been 
modified by chopper stabilization., The drift problem v/as 
negligible, for after initial alignment, checks on the 
less compllcatea units showed very little arift when 
noted daily. However, the more compllcatea systems, 
such as the first order lead/lag networks, and the sec- 
ond order generating system, which involves eight ampli- 
fiers in the units, required several checks in the course 
of solution. 

As will be seen in the attendant circuit diagrams, 
the operational amplifier potentiometers, which regulate 
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the magnitude of multiplication factors, are of the 
proportional type, ioe», one portion in the feedback 
loop of the DC amplifier, and the other is in the feed 
forwara loop® -^hese are normally set so that clockv»ise 
rotation of the potentiometer dial puts more resistance 
in the feedback, while reducing that in the forward loop, 
thus, causing multiplicationo In the K4DY second order 
however, the dial A 2 is wound oppositely, therefore div- 
ision is performed for clockwise dial rotationo This is 
necessitated by the form of the standard second order 
differential equationo 

Electronic clamping of the Integrating boxes back 
to zero potential is necessary, if the solution to the 
problem is not reached prior to the characteristic time 
of the integration units that are involved in an open 
loop solutiono If these units are connected in c losed 
loop, clamping is generally not requiredi however, due 
to d rifting DC levels, it is best to clamp all integrat- 
ing boxes, as a matter of courseo 

Calibration of the dial settings on the multipli- 
cation and DY boxes was difficult, since normal bridge 
resistance measurements would not give satisfactory re- 
sultso Rather, known voltages were fed into the boxes 
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being calibrated, and the dial set to give the desired 
value of output voltage „ This is not too time consuming, 
since other calibration methods on components require as 
much care and time, for accurate results » 

The feature of high speed solution and simultaneous 
viewing is the greatest advantage of this type of computero 
It is reaaily useful in simulating any control system, due 
to the many functions that can be generated, by a series of 

4 

interconnected black boxes, and a quick solution, or quick 
evidence of an improper solution, becomes apparent almost 
immediatelyo Since so many solutions, to physical systems, 
are in the form of the standard second order differential 
equation, then the K4DY box inherently is ideal for a quick 
approximate solutiono 

Method 

In setting up the pitch rate control system on 
the computer, there were some modifications required? 

In that the time constant of the unit lead box is 
set from zero to one, and it is oesired that all Tq he 
set at o 4ms (for ease in scaling) then the uncompensated 
open loop system equation of the form; 

p (2nd order equation) 
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is rewritten as: 



% 4_l)q„i 

P 

(2nd order equation) 

where 1 can be simulated by an integration boxo 
P 

The same type of problem arises from the compen- 
sated closed loop syscem form: 

q^ K(pTl+l) (pTp-fpq^^ where Ti)l 
(2nd order equation) . / 

and T 2 \ 1 

rewrite as: 

qo-K(p^l ) (pT2 +'l)T]^q^ 

Tl 

(2nd order equation) 

and KT-j_ becomes a new K’ , and p can be approximated by 
a derivative box» Phis derivative box, since a perfect 
derivative is almost impossible to generate, causes a 
time lag error of 1,5 seconds (problem time) in the paral- 
lel netv/ork in which it is operating. This leaus to an 
undesirable time lag in the computer simulation, xiealizing, 
however, that this particular shortcoming exists; permits 
its removal in the final analysis, Phis unwanted lag is 
not apparent in the unit lead box. 

The compensated system, has the open loop form: 
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qo=K(pT-l-l) (p + Z)qjji 

(2nd order equation) 

Vtihere 1 approaches one of the roots of the 2nd order 
T 

equation, and Z approaches the other rooto There is high 
gain in the forward loop, alsOc Then, to a close approxi- 
mation, the equation may be written in the final closed 
loop form as; 

qo^K- (pTj^+X)q„ 

; K' includes K and the T's 

(pT2 -f-1) that are determined^, as are 

and T 2 > according to the 
problem involvedo 

Thus, the equation can be mechanized in the closed 
loop by use of a multiplication box in series v/ith a unit 
lag and a unit lead box, with no feedback connection re- 
quired. 

Scaling 

Magnitude scaling was carried out, assuming that 
the peak signal voltage available from the square and sine 
wave generator was 13.5 volts and noting that 50v was the 
limiting factor of the DC amplifiers, and that the maximum 
input signal to the pitch rate system was from .1 to .3 
radians/second. For scaling purposes, the quantities E 
and X are introduced to the problem as intermediate quanti- 
ties. 
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E, the error. Is the difference between and 
X is the input to the unity gain second order systenio 

Fbr simplified presentation, the output gain of 
the system was changed to show the final results, in the 
photographs, in the proper relative aspect. That is, the 
output is scaled relative to the input. I’hus, if the in- 
put is 10 squares, and the output is 8 squares in steady 
state, then the final value is 80^ of the input, and no 
extra magnitude scaling is required of the reader. 

i'ime scaling for this problem was set up to be a 
special case, in that all units involving time were set 
for To equal to o4ms. Thus, 1 second of "real time" on 
the computer; in computer (problem) time is 2500 seconds 

A dual beam oscilloscope was used for final pre- 
sentation, in order that both input and output may be 
viev/ed simultaneously. 

Input Signals 

hesponse to a step function, to simulate violent 
maneuvers; and a ramp function, to simulate more gradual 
maneuvers, was used. i’he frequency of the square wave 
step function was 30 cps. In problem time, this is 42 
seconds/step, Ihe slope of the negative input ramp is 
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o 12 units of input/computer second. Each small square 
on the grid of the photographs represents lo67 seconds 
of problem time. 

Frequency Response 

The scope of this paper does not include experi- 
mental frequency analysis. However, in preliminary study 
of the problem it was found that the more complicated 
boxes have their own frequency characteristics. It was 
also noted that to make proper response for lead, or dif- 
ferentiating systems, the following must be done. Connect 
a system to s imulate a lead lag transfer function made 
up of Integrators, potentiometers and adders. Make the 
denominator term of the network far enough out on the 
real axis, so as to essentially leave a pure lead term 
in the final result, that gives the best frequency res- 
ponse. The unit lag boxes exhibited better frequency char 
acteristlcs. The second order box was accurate for the 
frequencies tested, and should adapt well to measurement 
of simulated sinusoidal gust type responses. 

Circuit Representation 

A more detailed study of the unit boxes, their 
functions and capabilities, along with a shorthand method 
of representation, is presented in later sections for the 
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following reasons: 

1. Since most present-day manufacturers of analog 
computers have all but abandoned the mooular technique; 
the internal circuitry is shown for comparison with present 
simplified ’’plug in” R and C component methodso 

2o To demonstrate, in block diagram form, the com- 
puter final circuit connections; to simulate the pitch 
rate system, model, and associated compensation devices.. 
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APPENDIX B-2 



Circuit and Block Symbols of the 
GAP/R Components o 



56 



/ 



[ 



-) 




i 




K5-A Adding Component 



Equation: y =■ -h xg + . . .X 




DESCRIPTION: The "Adder" computes the instantaneous sum 

of four or fewer input signals. By interconnecting Adders 
in tandem, any number of signals may be combined in sums 
and differences. V»ith N such Components, 3N -t-1 signals 
may be added. Each Adder features a calibrated dial where- 
by a steady voltage may be added in, amounting at the 
extremes to plus or minus 10 volts, or 20 % of the maximum 
signal excursion. This is convenient for manual addition 
of a constant term in an equation. 

Signals receive neither gain nor attenuation through 
the Adding Component. 
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K3-C Coefficient Component 




Block Symbol: 




DESCRIPTION: An input signal to this Component is 

automatically multiplied by a numerical factor, then 
carried directly and inversely to the two outputs. The 
factor may be set on a specially calibrated nonlinear 
dial, extending from 0 to 100. Technically, gains from 
a precise zero to infinity are within range, but, the 
best linear range is 0.1 to 10. 
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K3-J Integrating Component 



Equation ; *»’ y = 1 J xdt = x 

ToP 



Isec 




Block Symbol: 




DESCRIPTION: This Component, the "Integrator", computes 

an integral with respect to time of the voltage supplied 
to the lower input. The time-factor is 0.4 millisecond. 
For normal computing, in closed loops, the right-hand 
switch position, is chosen for maximum sensitivity; the 
left-hand switch position introduces a stabilizing net- 
work and provides for open-cycle operating when desired, 
as for the purpose of display alone. 

The integrated outputs are artificially and auto- 
matically returned to zero by means for a "clamping" sig- 
nal supplied to the upper input. 
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K3-D Differentiating Component 



Eauation; y a Tr^dx - T/^px 

dt 




Block Symbol: 



-f 

y 

-> — 




DESCRIPTION: Operationally, the "differentiator" 

performs Inversely to the Integrator. It computes 
time derivative of the input signal. The time -factor 
To is 0.4 milll-second. With the switch at the left- 
hand setting, as in this case, a lag is introduced 
having a time constant of about 50 microseconds; 
about l/lO the fixed time-factor of the Component. 

This setting is an approximate derivative, which 

did affect accurate representation showing as spurious 

lag. 
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Ko-E Augmenting Differentiator 



Equation: ± ysx -t- Tdx = (1+Tp)x 

dt 




y 




DESCRIPnON; This Component is related in a simple 
way to the K3-D Differentiator, since it passes the 
input signal and adds to it the first time derivative 
of that signal. The time-factor T, the derivative 
sensitivity, is continuously adjustable from zero up 
to 0.4 milli-second on a linear percentage scale. 

The operator embodied in inverse to that of the K3-L 
Unit-lag Component. 

Pbr certain purposes this Component may replace 
a combination of K3-A, K3-D, and K3-C. It represents 
the introduction of a derivative response. 
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K5-L Unit-Lag Component 



Equation: y ♦ s(l + Tp)y=+x 

dt 




Block Symbol: 



X > 




DESCRIPTION: The "Unit-lag” performs a first-order 

lagging operation. By means of its 0-100 dial, the 
time constant T may be adjusted from zero up to a 
maximum of 0.4 millisecond. The final, or long-term, 
sensitivity is unity. This type of operation or 
dynamic characteristic may be obtained by a simple 
loop involving a K3-A, a K3-J, and a K3-C, which 
loop also affords several other characteristics. 
However, the K3-L Component provides a compact and 
convenient version of this frequently recurring 
operator. 



62 



K4DY Dynamic Component 



Equation: AgT§ ^ ^ A^To ^ AqY = x(t) 

dt^ dt 




* 7-7 



DESCRIPTION; This box combines the Tunctions of eight 
assorted A,G, and J boxes into one simplified unit. The 
standard second order equation satisfied by this unit, 
is an intrinsic part of most actual physical systems. 

Tq is equal to .4 milliseconds, and Aq,A 2 »Ao are directly 
calibrated on 0-100 dials of the Coefficient type. Dial 
A 2 is wound inversely, as previously described. 
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APPENDIX B-3 



GAP/R Block Connection Diagrams and Dial 



Settings, 



for the Pitch Rate System Simulation, 
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FLIGHT CQNDITICM A 




FLIGHT CONDITION B.C 




FLIGHT COMDITICN D.E.F.G 




Computer Block Diagrams for the Compensated Pitch Rate 
System, During the Various Flight Conditions. 
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FLIGHT CQMDITIOM A,C 




to scope 



FLIGHT CQNuITiOM 




to scope 

Computer Block Diagrams for Uncompensated Unity Feedback 
during the Various Flight Conditions, Pitch Rate System. 
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Wji s 3*5 radians/second 

r r .7 

Dial Settings for the Model 

^1>C2-1 Cg»4,9 04 , 055 * 12.25 Ji, Jg = *4rasec . 



Ooraputer Diagram and Dial Settings for the Model 



67 










I 



Unit Dial Settings for Compensated System 



Unit 

PlightN. 

C ond i ti on^\^^ C ^ 


C 2 C3 


E 


CM 

< 


Al 


Aq 


A 


,185 


.500 1 


.125 


ol 


2.8 


4 


B 


o936 




.126 


,123 






C 


.999 




.125 


,125 






D 


J — 1 
0 

00 


4.37 


.124 


5 


4.56 


1.03 


E 


8.43 


4.82 


.125 


5 


4.43 


1.03 


P 


00 

0 

0 

CO 


2.32 


,125 


5 


4.71 


1.44 


G 


7.98 


6.15 


.125 


5 


4.73 


1.43 



a» K4DY and K3D set for To — o0004 sec, 

bo A 2 , and Aq are the dial settings of the 
K4DY 2nd order equation module® 



68 




I 





Unit Dial Settings for the Uncompensated 
Systemc (with unity feedback) 



X^nit 

Fllgh^X 

Condition^ 




^2 


^3 


C4 


E 


^2 


Al 


Ao 


A 


o925 


o925 


o69 




o935 


2o5 


5o25 


lo21 


B 


loll 


loOO 


o073 


o094 




o500 


o300 


o092 


C 


o925 


o925 


ol67 




o325 


2o5 


15ol 


2o32 


D 


loll 


loOO 


o018 
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a» K4DY and K3J are set for Tq = <,0004 sec„ 

b„ Agj A]_, and A^ are the dial settings of the 
K4DY 2nd order equation module <> 
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APPENDIX C 



Transient rieponse of the Pitch Rate Systemo 
1. General information and definitionso 
2o Step input, simulating violent raaneuvero 
3o Ramp input, simulating gradual raaneuver. 
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GENERAL INFORMATION AND DEFINITIONS 



lo Theoretical step response of the compensated pitch 
rate system, determined by inverse transform methods o 



FLipjht Condition 
A 
B 
C 
D 
E 
F 
G 



Step Response 
o 965 + ,0222 e “^‘’^t 
qo = .945 + o0245 e “8»155t 

qQ=o990- .00497 e -8“04t 

- ADRI- -8o5t 

qo^ .768+ .092 e “°^^^+.0545e 

qo= .84 +.0357 e +.073 

q^^ ,6 + .268e‘“^25t4.^0552e“®’^^ 

qo= .761+ .0587e .067e“®° 



2, The uncompensated, unity feedback, system reduces to 

a- closed loop type 1 servo system. The final value theorem 
predicates that the final amplitude ratio of output to input 
is equal to 1. 

3. Each small block on the photographs is equal to 1.67 
seconds of problem time. 



4. For the ramp input function; the output ramp is ex“ 
trapolatea, ana the ''aynamic error" is estimated. 
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Dynamic Error (DE) is defined as the difference 
between the output response level at t, and the forcing 
function point for the same instant: 
or in equation form: 

^^E)t]^= <l(out) ti “■ ^(out) qs ti 
where qs = quasi-static, the estimate'of the output 

level rate, for quasi-static conditions • 

(DE) is then nondimensionalized as ^^/iMPUT 

5o For the step input function, the "response time" is 
determinedo 

Response time (RT) is defined to be the time required 
for the output response to change from an initial arbitrary 
initial dynamic error to within five percent of the forced 
output response, or to within 95% of its steady state valueo 
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(a) Compensated 
RT ,835 seconds 
DE/Input .05 



(a) Uncompensated 
RT 5.85 seconds 
DE/Input .1 
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(B) Compensated 
RT ,167 seconds 
DE/Input .1 




(B) Uncompensated 
RT 10.85 seconds 
DE/Input .581 
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■ (C) Uncompensated 

RT 10.85 seconds 
DE/Input ,561 
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I 

I 



I 




1 



■ (D) Compensated 

RT 5.0 seconds 
DE/Input .474 

■ (D) Uncompensated 

RT 13.4 seconds 
DE/Input .5 
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■ (E) Compensated 

RT 7.51 seconds 
DE/Input .5 

■ (E) Uncompensated 

RT 10.85 seconds 
DE/Input .208 



77 







(F) Compensated 
RT 4.58 seconds 
DE/Input .465 



■ (F) Uncompensated 

RT 15.08 seconds 
DE/Input ,615 
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■ (G) Compensated 

RT 6.51 seconds 
DE/Input .85 

■ (G) Uncompensated 

RT 17.5 seconds 
DE/Input .80 
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APPENDIX D 



Transient and Frequency Response of 
the Combined Model and Pitch Rate 
System* 



80 




1 



I 



[ 




The perrormance function of the model in 
closed loop is: 

q in 

p2f 2fWnP+wg 

vi/here: Wn= 3o5 radians/second 

f = .1 
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Combined System 



RT 1.1 seconds 



DE/input 



0 
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